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COMPUTER  GRAPHICS:  TWO-  AND  THREE-DIMENSIONAL  CLIPPING 


L INIHOOUCnON 

cupping  is  a procedure  used  to  determine  the  proper  line  segment(s)  of  a given  Une,  such  that  the 
resulting  segment(s)  contain  nq  points  exterior  to  a {dane  that  serves  to  define  a bound^.  Clipping  is 
em{doyed  to  limit  the  amount  of  work  a disiday  device  must  perform  by  drawing  oidythe  visible  parts  of  a 
Une.  This  is  done  by  uting  a dipping  algorithm  to  compute  only  the  line  s^ments  interior  (or  exu^r)  to  a 
defined  region  that  has  been  selected  for  disiday.  The  resulting  list  of  lines  represents  a substantiaUy  siiudler 
amount  of  work  for  the  disjday  hardware  to  process  in  a given  period  of  time. 

A good  dipping  algorithm  quickly  rqects  Unes  that  Ue  outside  the  viewing  area,  and  since  speed  is 
always  desirable,  algorithms  that  lend  themselves  to  parallel  (or  pipeUne  imfdementatian),  in  either 
hardware  or  firmware,  are  sou^t. 

Let  us  initiaUy  limit  our  discussicm  to  two-dimendonal  dipping.  A later  section  of  tius  paper  wfl 
exidore  three-dimensiotud  dipping. 


n.  TWO-DIMENSIONAL  CUPPING 

Typically,  disiday  regions  have  been  diosen  as  being  rectangular  (Rgure  1).  The  actual  number  of 
sides  is  not  of  importance  to  the  development  of  a dipping  algorithm,  alftou^  it  does  affect  performance 
when  finaUy  implemented,  dewing  windows  are  always  diosen  however,  as  convex  sets  in  a lUbert  Space. 


^R 


The  standard  properties  of  lUlbert  spaces,  as  well  as  the  set  concepts  referred  to  in  this  paper,  are 
found  in  most  texts  on  linear  analysis. 

Making  use  of  these  properties,  the  foUowing  is  dear:  A straight  Une,  intersecting  the  interior  of  a a 
convex  set,  can  intersect  the  boundary  of  the  set  in,  at  most,  two  places.  In  the  event  the  ponvex  set 'is 
dosed  and  bounded,  the  straight  line  wffl  intersect  it  in  exactly  two  places. 

At  this  point,  we  introduce  the  idea  of  ‘‘inward”  normals.  That  is,  we  say  n is  an  inward  normal  for 
any  boundary  point  X of  S if,  for  any  other  point  Y in  S, 

n * (Y-X)>  0.  where  denotes  the  inner  product  operation; 

With  this  representation,  we  have  defined  a natural  geometry  ovpr  the  set  S. 

For  the  set  S only  four  unique  inward  normak  exist.  AU  othen  are  mathematically  equivalent  to 
these  four. 

Our  particular  dipping  problem,  tiut  of  determining  the  points  of  intersection  between  a particular 
Une  segment  and  the  set  S can  foen  be  partidly  resolved  by  examining  tiie  parameterized  behavior  of  the 


W, 


Figure  I.  Rectangular  dbplty  region. 


line  lepnent  with  the  four  inward  normals  to  S.  In  particular,  if  f it  a boundary  point  of  S for  whidt  n it  an 
inward  sector  and  P(T)  is  a paiameteriiation  of  the  line  segment  in  question,  then,  for  particular  values  of 
T. 

n-  [P(T)-fl>0 

ipippi^  P(D  is  “pointed"  toward  the  interior  of  S with  respect  to  at  least  one  of  the  four  notmalt, 
n-lP(T)-fl  *0 

implies  P(T)  is  “pointed”  parallel  to  the  0ane  containing  f and  perpendicular  to  n (that  ta,  parallel  to  this 
particular  bound^  of  S)  and 

n-  [P(T)-fl<0 

iitqdies  P(T)  is  “pointed"  away  from  the  set  S. 

So  far  we  have  mentioned  only  that  any  line  piercing  S does  so  at  exactly  two  points.  It  is  also  dear 
that  these  two  points  of  intersection  do  not  lie  on  the  same  target  fdanes  to  the  boundary  of  S (a  plane  not 
completely  containing  a line  may  be  intersected  by  it  at  most  once).  These  facts  give  rise  to  our  first 
irrrportant  result: 

if  f is  a p<*rt  in  a particular  pdane  bounding  S for  vdiich  n is  an  inward  normal,  then  the  scalar 
equation 

n‘  [F(T)-fl  =0 
has,  at  most,  one  solution. 

Thus,  there  arc  exactly  two  possibilities;  first,  P(T)  U parallel  to  some  pdane  lying  along  die  boundary 
of  S.  Second  P(T)  pierces  that  plane  for  some  value  of  T.  Now,  suppose  we  want  to  talk  about  drawing 
only  a line  segment: 

:P(T)  = Po+(Pi-Po)TTe(0,l] 
h h the  restriction  of  T to  [O,!]  that  makes : P(T)  a line  segment. 

With  regard  to  the  set  S,  either: 

1.  Ihe  line  segment  is  wholly  contained  in  S,  or 

2.  The  line  segment  is  partly  contained  in  S,  or 

3.  The  line  segment  is  no  part  of  S at  ail. 

Before  introducing  the  prrimary  result,  we  need  the  foDowing  definitioru: 

Define  Sj,  as  fire  set  of  all  planes  bounding  S,  that  Is,  intersecting  only  die  boundary  of  S. 

For  each  (,eS(,,  define  n(b)  u the  set  of  positive  normals  drawn  from  b with  respect  to  S. 

The  Primary  Result 

The  set  of  an  T sudi  that: 

1.  Te(0,l) 

2.  n-  lP(T)-fI>0 

Where  for  every  beSj,,  n e n(b),  and  f b an  arbitrary  vector  in  b,  completely  determines  die  line  segment 
contalnineot  far  S. 

Careftd  exandnation  of  the  flowchart  (Rgure  2)  reveals  the  underlying  geometry  involved  in  this 
technique.  The  tanpdementatkm  lends  itself  to  swift  line  sepnent  rejection. 
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The  previous  result  is  great  for  theory,  but  the  practice  is  somewhat  more  comi^ex. 

1 . How  many  planes  are  involved  in  the  convex  object?  If  the  result  were  large,  our  procedure  would 
be  tedious. 

2.  Are  there  any  ways  to  “set  up”  the  proUem  to  facilitate  the  sdution? 

3.  Are  there  any  timesaving  ways  to  accept  or  reject  segments? 

It  should  be  dear  that  there  is  no  one  fastest  way  to  examine  any  vector  segment  under  all 
circumstances.  For  example  (Figure  3),  examine  that  rectangle  again: 


AUTOMATIC 

REJECTION 


Figures.  Probabilistic two^limensional dipping. 


V'  abe  an  ellipse  and  circumscribe  an  ellipse  and  test  P[T)  for  either  immediate  acceptance 

or  reje  ju.  oi  the  event  the  test  fails,  the  area  we  must  further  test  is  considerably  smaller.  This  kind  of 
test  procedure  is  one  type  of  conditional  test,  where  we  maximize  the  probability  of  an  accurate  selection. 
Now,  we  consider  the  general  formulation  of  the  orientation  vector  technique. 

Let 

P(T)  = Po+[P,  -Po]TTe[0,l] 

P(T)  can  be  made  a line  by  allowing  Te(— <»,<»). 


Figure  4.  Rectangular  dtaplay  region  conventioiis. 
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Hie  flve  conditions  become 

lP(T)-fi]-ni>0  i = 1^2,3,4 
orPo  • ns+T[P, -Pol  ■ ni>0  i*l,2,3,4 
orTfP, -Po]  nj>[fi-Pol-ni  i = l,2,3,4 
tndT6[0,l] 
letting 

Pi  -Po  = D (Die  called  the  Directorix  of  the  vector  P(T)) 
and  Wj  *Po-fj(Wj  is  the  ithjAane  “weight”  vector),  i = 1,2,3, 4 
Ihen  we  can  write: 

TOni  + Wi*ni>0  i = 1,2, 3,4 
Now  we  have  our  algorithm; 

For  a K-sided  convex  object  S with  inside  normals: 

n, , nj, ...,  njj 

Points  on  each  {dane: 

fi, fa  fjj 

Let  P(T)  be  the  parameterization  of  the  line  running  from  Po  to  Pi 
Define: 

Wi  = Po-fi,D  = P, -Po 

Then  the  K + 1 conditions  for  any  portion  of  the  luie  segment  being  in  S are: 

TD-iii  + Wj-nj>0  i=l,2 K 

Te[0, 1] 

Let  the  set  of  such  that  each  ith  case  above  is  satisfied: 

Tj=  |t  |TD-nj  + Wj-ni>0| 

Then  for  eadi  Tj,  Tjc  (-“,«») 

LetT=MTj 
i = l 

1henT=  |TlP(T)eS^ 

Rnally,T*  = Tn[0,ll 

Using  the  Tektronix  4051  Graphic  Computing  System,  an  implementation  of  this  algorithm  was 
programmed  (Rgure  5).  The  program  and  sample  results  follows  (Figure  6): 

In  this  program: 

Pbx  = Xl  Plx  = X2 

POy*Yl  Ply=Y2 

Figures  7 and  8 show  examples  of  segments  drawn  and  rejected,  respectively. 

We  can  introduce  efficiency  into  the  computation  by  first  evaluating  directly  for  the  values  of  T. 
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100  INIT 
110  PAGE 

120  WINDOW  0,22,0,22 
130  VIEWPORT  65,130,0,65 
140  DIM  F(4,2)J4(4,2) 

1 50  REM Define  window 

160  Wl=7 
170  W2=18 
180  W3=10 
190  W4=17 

200  REM Draw  axis 

210  MOVE  0,0 
220  DRAW  22,0 
230  MOVE  0,0 
240  DRAW  0,22 

250  REM Draw  Window  boundary 

260  MOVEWl,W3 
270  DRAWW2,W3 
280  DRAWW2,W4 
290  DRAWW1,W4 
300  DRAWW1,W3 

310  REM Enter  line  definition 

320  PRINT  @32,21,0,90 

330  PRINT“Enter  xl  ,y  1 , x2,y2  : ” ; 

340  RWT  XI  ,Y1,X2,Y2 

350  REM Nomuds  defined 

360  N(1,1)=W2-W1 
370  N(l,2)=0 
■380  N(3,l)=-N(l,l) 

390  N(3,2)=-N(l,2) 

400  N(2,l)=0 
4.10  N(2,2)=W4^W3 
420  N(4,1)=-N{2,1) 

430  N(4,2)=-N(2,2) 

440  REM Window  boundary 

450  F(1,1>=W1 
460  F(1,2)=W3 
47''  F(2,1)=F(1,1) 

480  F(2,2)-F(l,2) 


490  F(3,1)=W2 
500  F(3,2)=W4 
510  F(4,1)=F(3,1) 

520  F(4,2)=F(3,2) 

530  REM Algorithm 

540  D1=X2-X1 
550  D2=Y2-YJ 
560  L=0 
570  U=1 

580  FORM  TO 4 
590  W1=X1-F(I,1) 

600  W2=Y1-F(I,2) 

610  D=D1*N(I,1)+D2*N(I,2) 

630  W=W1*N(I,1)+W2*N(I,2) 

630  IF  D=0  THEN  690 
640  T=-y//D 
650  IF  D>0  THEN  710 
660  IF  T<0  THEN  860 
670  U=TMINU 
680  GO  TO  730 
690  IF  W<0  THEN  680 
700  GO  TO  730 
710  IF  T>1  THEN  860 
720  L=TMAXL 
730  NEXTl 

740  REM Segment  drawn 

750  IF  L=>U  THEN  860 

760  X4=X14(X2-X1)*I 

770  Y4=Y1+(Y2-Y1)*I 

780  X5=X1+<X2-X1)*U 

790  YS»Y14(Y2-Y1)*U 

800  MOVEX4,Y4 

810  DRAWX5,Y5 

820'  PRINT  @32,21:0,30 

830  PRINT  “SEGMENT  DRAWN” 

840  END 

850  REM Segment  rejected 

860  PRINT  @32,21:0,30 

870  PRINT  “SEGMENT  REJECTED” 

880  END 


Figure  5.  Sample  computer  program. 
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(1)  T(P,^-Po^)(Wr-Wl)>(Wl-Po^)(Wr-Wl) 

(2)  T (P, y - Poy)  (W^Wg)  > (Wg-Poy)  (Wt-Wb) 

(3)  T (P.^  -P„^)  (Wl-Wr)  > (Wr-P„^)  (Wl-Wr) 

(4)  T(P,y-Poy)(WR-WT»(Wp-Poy)(WB-WT) 

(5)  Te[0,l] 


Assiu^g  that  KT)  is  not  parallel  tn  any  side  of  the  rectan^e,  then  we  can  sdve  for  the  values  of  T 
so  that  each  equation  (1  through  4)  is  satisfied  identically. 


m.  THREE-DIMENSIONAL  CUPPING 

^ce  the  algorithm  was  not  restricted  to  two  dimensions,  the  basic  result  applies  in  the 
Aree-dimensional  case.  We  will  examme  a particulai  case.  The  set  of  points  chosen  for  perspective  purposes 
is  centered  with  the  axis  center  at  the  eye  point,  looking  along  the  2 axis.  The  set  of  points  is  given  by  the 
set  product 

|Wl.Wr|  X |Wb,Wt|  X |WH,WYf 

The  region  in  space  R,  now  being  viewed  as  our  convex  set,  is  given  by  the  product 

[Wl,Wr]  X [Wb,Wt1  X [Wh.  Wy] 

In  the  2-D  case,  there  were  four  equations  in  T and  one  additional  constraint.  In  this  3-D  case,  there 
are  six  equations  in  T and  that  same  constraint. 

There  are  six  bounding  planes,  four  of  adiich  contain  the  point  [0,0,0].  The  other  two  are  the  hither 
and  yon  planes,  which  are  handled  as  special  cases  (Figure  9).  Symmetry  conditions  imply  the  foUowiig: 

^L=-^R  Wb*-Wx 


V,=[Wr.Wt.Wh] 

V,  =I-Wr.+Wt,Wh1 

V3=.[-Wr*-Wt.Wh1 

V4  = [+Wr,-Wt.Wh1  Figure  9.  ThreeHUmensioiial  viewing  window. 

These  four  vectors  are  the  vectors  along  the  intersection  of  the  four  bounding  planes  to  the  first 
perspective  plane. 

n,  = VI  X V,  = [0,  -2WrWj,,  2WrWj] 
n,  = V,  X V3  = PW^Wh.  0,  2WrWt1 
n,  = V3  X V4  - [0.  2WrWjj,  2WrWt1 
n«  = X V,  = [-2WtWh.  0,  2WrWt] 

n,  = (Hither)  = [0,0, 1 ] 
n«  = (Yon)  = [0, 0,-1] 


D=  P.  -P, 

Wl  = -f,  + Po 

W.  - +P„ 


W,  = + P, 


W,  = +P, 


w.  =+P„ 


The  same  rules  apply  to  O'nj  and  Wj ' nj  as  before. 

The  primary  formula  is  T = ~*^i  **H  for  each  i,  giving 

D • nj 


WT-Poz+W^Poy 

“ O'.,  % - (P,y  -Poy)WH 


WhP,,+WrPo, 


0*.x  + -Po,)Wr 


z ”z 


WnPoy  + WTPo, 


(P.„  - PoJWh  + (P.^  - Po,)Wt 


• y »y 


T* 


-WrPoz  + WhPox 


(P.,  - Po,)Wr  - (P.,  - Po,)Wh 


T, 


-Po, 


T. 


^ -"«z 


These  values  of  T,  define  the  intersection  points  of  P(T)  with  each  bounding  plane. 

It  can  be  teen  from  this  analysis  that  while  the  computation  of  T for  the  3-D  case  may  appear  to  be  a 
lengdiy  computation,  each  value  of  T can  be  computed  in  a separate  parallel  process. 

Our  discussion  to  far  hat  been  limited  to  deteimiidng  a line  segment  that  is  either  rejected  or  dipped 
to  fit  inside  a convex  set  S.  Our  output  (T^.Ty)  can  also  be  used  to  draw  the  complement,  that  is  the  line 
segment(t)  exterior  to  a given  convex  set  S. 

The  poMibilitiet  ate  summarized  in  Rguie  10. 
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Hole 

Exterior  Dipping 

Hole. 

Interior  Gippinf 

Normal  exit  draw  P(Tj^),P(Tjj) 

draw  P(Ti;:P(Ty) 

Tl=^),Tu=1  omit 

TfTu 

’'l'Tu 

omit 

X 

Abnormal  exit 

Omit 

Tl=o.Tu  = i 

Figure  10.  Window -hole  dipping  summary- 


This  introduces  the  possibility  that  an  Algol,  PL/ 1 or  Pascal  implementation  of  this  algorithm  might 
be  implemented  to  first  apply  exterior  dipping  to  a ‘Hole’  contained  within  the  convex  window. 

If  we  define  a function  which  performs  P(T)  ca)^  'EG,  the  flow  chart  (Figure  2)  previously  drawn 
can  be  extended  as  in  Figure  1 1 . The  output  of  SEv.  . scaled  appropriately  for  a display  device  using 
integer  coordinates. 
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Exterior  Gip 


Interior  Gip 


Normal  Exit 


Abnormal  Exit 


exterior 

clip 


exterior 


CALL  SEC 
(L.1X2.1Y2) 
CALL  SEC 
(U.IX3.1Y3) 


CALL  SEC 
(0,0,1X1.1Y1) 
CALL  SEC 
(1.0, 1X4, 1 Y4) 


CALL  SEC 
(O.O.IXLIYl) 
CALLS^ 
(L,IX2,1Y2) 


CALL  SEC 
(U,IX3.IY3) 
CALL  SEC 
(1.0,IX4,1Y4) 


Figure  II.  Extended  dipping  flowchart, 


WINDOW  Wl  W^  Wg  Wy  Wj^  Wy  E 
HOLE  Hl  Hg  Hg  Hy  Hy 

VIEWPORT  Vl  Vg  Vg  Vj 


Where:  Wl  < Hl  Wg  > Hg 
Wg  ^ Hg  Wy  ^ Hy 
Wh  < H|^  Wy  > Hv 


Figure  12.  Wiiulow-hole  dipping  perspective. 


IV.  CONCLUSION 

Implementation  of  true  three  dimensional  interior  and  exterior  clipping  gives  a disi^ay  processor  great 
flexibility  in  creating  a troly  useful  representation  of  complex  shapes  short  of  hidden  Hnc  elimination.  By 
careful  placement  of  the  boundary  planes,  the  viewer  will  be  able  to  sectirm  shapes,  remove  obscuring 
surfaces,  or  look  through  one  surface  to  another  one  in  the  interior.  These  techniques  will  form  part  of  the 
necessary  tools  visual  displays  win  find  necessary  to  provide,  as  computer  graphics  seeks  to  aid  in  the 
visualization  of  objects  presently  being  dealt  with. 
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